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LINEAR  THEORY  OF  THE  E  X  B  INST  ABILITY  WITH  AN 
INHOMOGENEOUS  ELECTRIC  FIELD 

Z.  INTRODUCTION 

Ail  la port ant  Instability  associated  with  the  structuring  of 
Ionospheric  plaanas  (e.g.,  high  latitude  F  region  and  barlun  clouds)  Is 
the  £  x  \  Instability,  also  known  as  the  gradient  drift  Instability.  The 
Instability  is  an  Interchange  Instability  which  can  occur  in  an 
lnhoaogeneous ,  weakly  collislonal,  magnetized  plaaaa  that  contains  an 
aablent  electric  field  orthogonal  to  both  the  aablent  magnetic  field  and 
the  density  gradient.  A  simple  physical  picture  of  the  Instability 
mechanism  is  shown  in  Fig.  1.  We  consider  a  plasma  such 

A  A 

that  B  -  B  e8,  £  «  E  ey,  n  -  n(x)  with  3n/3x  >  0  and  ven/ae  <<  vin^ai 
«  1  where  is  the  collision  frequency  between  spedee  a  end  neutrals, 
end  Qa  is  the  cyclotron  frequency  of  species  a.  Upon  this  plasms  we  Impose 
a  density  perturbation  dn  ~  6n  sln(k^y)  as  shown  In  Fig.  1.  The 
influence  of  &  on  the  plasma  is  to  cause  (1)  the  electrons  and  ions 
to  |  *  l  drift  in  the  z  direction  end  (2)  an  Ion  federsen  drift  In  the  y 
direction.  The  latter  effect  Induces  a  space  charge  perturbation  electric 
field  denoted  by  6£.  The  response  of  the  plasma  to  this  perturbed 
electric  field  Is  to  drift  with  a  velocity  -  c6E  x  g/B2.  For  the 
configuration  shown  in  Fig.  1,  6^  causes  the  "heavy"  fluid  perturbation  to 
fall  Into  the  "light"  fluid  (region  I),  and  the  "light”  fluid  perturbation 
to  rise  into  the  "heavy"  fluid  (region  II)  -  the  classic  Interchange 
phenomenon.  Of  course,  if  the  direction  of  3n/3 be  or  were  reversed  then 
the  density  perturbation  would  be  damped. 

The  original  study  of  the  Ixl  Instability  urns  by  Simon  (1963)  and 
B»h  (1963),  who  applied  it  to  laboratory  gas  discharge  experiments. 
Subsequent  to  these  first  Investigations,  a  considerable  amount  of 


nmreh  bw  been  dwotad  to  npUinlsi  ionospheric  phenoaena  based  upon 
this  instability  (Union  and  Wortaan,  1970  and  references  therein;  Siaon, 
1970;  Polk  end  Bserendel ,  1971;  Perkins  et  al.,  1973;  Zabueky  at  al., 
1973;  Shiau  and  Siaon,  1972;  Ptrkins  and  Doles,  1975;  Seanaapleeo  et  al., 
1979;  Chsturvedl  and  Ossakov,  1979;  Keaklnen  and  Osaakow,  1982).  Two 
areas  of  present  interest  concerning  the  instability  are  barion  cloud 
etrlstlotts  (see  for  essnple  the  review  papers  Ossahow  (1979)  and  Ossekow 
et  si.  (1982),  sad  the  refsrences  therein)  and  the  structuring  of  plaana 
“blobs"  in  the  high  latitude  F  region  (Vickrey  et  al.,  1980;  Btsklnen  end 
Oseahow,  1982). 

Is  purpose  of  this  paper  is  to  Resent  s  general  theory  of  the  {  x  | 
Instability  which  considers  an  aablent  electric  field  at  an  arbitrary 
angle  to  the  density  gradient,  and  allown  the  electric  field  component 
parallel  to  the  daaslty  gradient  to  be  lnhoaogeneous .  gone  sspecte  of  the 
prehlen  have  boon  treated  by  Parkins  et  al.  (1973)  and  Perkins  and  Boles 
(1975).  IWrkins  and  Bales  (1975)  node  the  In portent,  discovery  that  the 
sheared  velocity  floe  (resulting  trm  •*  lnhoaogeneous  electric  field 
parallel  to  the  density  gradient)  can  stabilise  the  instability. 
Fttrtheraore,  short  wavelength  nodes  are  preferentially  stablUsed  over 
longer  wavelength  nodes.  The  work  of  Psrkins  end  Boles  (1975)  considered 
the  strong  collision  Unit  (v^Q  »  a),  assuned  a  specific  density  profile 
aasnable  to  analytical  theory,  and  is  valid  only  in  the  short  wavelength 
regiaa,  l.e.,  kyL  »  1  where  L  is  the  scale  length  of  the  boundary 
layer.  The  present  study  extends  the  theory  of  Psrkins  and  Boles  (1975) 
by  renewing  these  restrictions.  Hasely,  we  derive  a  differential  situation 
which  describes  the  node  structure  of  the  I  x  |  instability,  ion  inertia 
effects  era  Included  so  that  the  ratio  v^/a  is  arbitrary.  Moreover,  we 


lolvt  this  equation  ganrlealljr  so  that  arbitrary  danalty  and  electric 
field  prof  ilea  can  be  considered,  and  the  reglaa  k^L  <  1  can  be 
Investigated  self-consistently. 

The  principal  results  of  this  work  are  the  following* 

1.  The  baalc  conclusions  of  Parkins  and  Doles  (1975)  are  verified 
ttuaarlcally.  Specifically,  the  marginal  stability  criterion  they  derive 
analytically  agrees  wall  with  our  numerical  result. 

2*  The  marginal  stability  criterion  la  weakly  dependent  upon  the 
magnitude  of  V- 

3.  The  stabilisation  mechanism  is  associated  with  the  x  dependent, 
Doppler-shifted  frequency  «  -  ky?y<x),  where  Vy(x)  ■  -  clx<x)/B,  and  not 
velocity  ehear  terms  proportional  to  9yh  ot  *2Vy/ht2. 

4.  When  i(s.)  >  I.,  where  s  10  the  pooitloe  about  Aich  the  mode 

la  localised,  the  moat  unstable  modes  have  kJL  <  1. 

7  m 

The  organisation  of  the  paper  Is  as  follows*  Da  Section  II  we  derive 
the  mode  structure  equation  for  the  (  *  }  Instability.  la  Section  III  we 
present  both  analytical  and  numerical  results  baaed  upon  this  aquation, 
finally,  in  Section  If  we  sunmarlsa  our  results  and  discuss  applications 
to  ionopherlc  phenomena,  l.e.,  barium  cloud  strlations  mad  high  latitude  P 
region  irregularities. 
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II.  THEORY 

The  equilibrium  configuration  used  in  the  analysis  is  shown  in  Fig. 
2.  The  aablent  Magnetic  and  electric  fields  are  in  the  a  direction  and 

A  A  A 

the  xy  plane,  respectively,  where  B  -  B  e_  and  E  »  E  (x)  e  +  E  e  .  The 

^  a#  x  x  y  y 

electric  field  in  the  y  direction  is  constant,  while  the  electric  field  in 
the  x  direction  la  allowed  to  be  a  function  of  x.  This  gives  rise  to  an 
inhoaogeneous  velocity  flow  in  the  y  direction,  l.e.,  Vy(x)  -  -cEJt(x)/B. 
Hie  density  is  taken  to  be  inhomogeneous  in  the  x  direction  (n  -  n(x»  and 
temperature  effects  are  Ignored. 

The  basic  assumptions  used  in  the  analysis  are  as  follows.  He  assume 
that  the  perturbed  quantities  vary  as  <p  ~  6p(x)  exp  [I(kyy  -  «t)J,  where 
ky  is  the  wave  number  along  y  direction  and  u  •  u»r  +  ly,  implying  growth 
for  Y  >  0.  The  ordering  In  the  frequencies  is  such  that  u  «  and 
vln«  &1  (the  F  region  approximation),  where  vln  is  the  ion-neutral 
collision  frequency  and  Qi  is  the  ion  gyrofrequency.  He  neglect  terms  of 
order  «/Qt  and  \n/\»  But  retain  terns  of  order  vin/».  He  Ignore  finite 
gyroradlus  effects  by  limiting  the  wavelength  domain  to  kr^  «  1,  where 
r^  is  the  mean  ion  Larnor  radius.  Hi  neglect  perturbations  along  the 
swgnetlc  field  (k(  ■  0)  so  that  only  the  two-dimensional  mode  structure  in 
the  xy  plane  is  obtained.  He  retain  ion  inertial  effects,  thereby 
including  the  ion  polarisation  drift,  but  Ignore  electron  inertia. 

A  key  feature  of  our  analysis  is  that  a  nonlocal  theory  is 
developed.  That  is,  the  mode  structure  of  the  potential  in  the  x 
direction,  the  direction  in  which  density  and  the  flow  velocity  are 
assumed  to  vary,  is  determined  by  a  differential  equation  rather  than  an 
algebraic  equation  obtained  by  Fburler  analysis.  This  is  crucial  to  the 


5 


analysis  since  Perkins  and  Doles  (1975)  have  shown  that  a  nonlocal 
analysis  is  necessary  to  demonstrate  the  stabilising  Influence  of  velocity 
shear,  due  to  the  inhomogeneous  electric  field  parallel  to  the  density 
gradient . 

The  fundamental  equations  used  in  the  analysis  are  continuity  and 
momentum  transfer: 


TT  *  7  ’  <"aV 


°*-f  a  +  U'» 


«■  +  i  h  * ») 


V.  E. 
in  n 


where  a  denotes  species  (e:  electrons,  i:  ions)  and  other  variables  have 
their  usual  meaning.  Rote  that  electron  Inertia  terms  are  neglected,  but 
ion  inertia  terms  are  Included,  and  that  the  neutral  wind  is  assumed  to  be 
sero.  The  equilibrium  drifts  are 


X*  -  -  cEx(x)/B  ey  (4) 

Xt  *  ((vin/0i)  cEx(x)/B  +  CVB) 

+  (-cEx(x)/B  +  (vln/fl1)  cEy/B)  ey  (5) 

where  we  have  chosen  a  reference  frame  such  that  Vx  ■  Vx  -  cEy/B, 

-  eB/m^c  and  E'(x)  •  3Ex/3x.  A  relationship  between  n(x)  and  ^(x)  can 
be  derived  by  assuming  V  •  £  «  V  •  [nCV^  -  X*)!  *  ®  which  reduces  to 


v  •  (n  -  0  (6) 

Equation  (6)  loads  to 

n(x)  ((Y^/Oj)  cEx(x)/B  +  (cE'(x)/BBi)  cEy/B)  -  constant  (7) 

where  we  take  the  constant  to  be  the  LHS  of  Eq .  (7)  evaluated  at 
x  ■  -  Thus,  by  specifying  the  density  profile,  the  electric  field 

profile  Eg(x)  can  be  determined  from  Eq.  (7).  Of  course,  if  there  are 

sources  and/or  sinks  in  the  plasma  such  that  V  •  (n  V^)  *  °»  then  Eq. 

(7)  Is  not  applicable. 

We  now  consider  a  linear  perturbation  analysis  of  Eqs.  (l)-(3).  We 

assune  ®0  "  n0  +  fina»  Xa  “  Xa  +  5Xa  and  E  -  £  -  V$  where  $  is  the  perturb¬ 

ed  electrostatic  potential.  Using  Eqs.  (2)  and  (3),  we  obtain 


6Vex  "  “  lky*  <C/B>  (8) 

SVey  *  (c/B>  <9> 

-  (-ikyd  -  v-y/a1)  *  +  i(«/n1)  -  (Vlx/a1)  ♦")  (c/b)  (io) 

6Viy  “  +  1  ♦  +  (1  “  ikyvlx/ni)  ♦')  <C/B>  (H) 


where  m  -  w  +  ivln  -  kyVly,  V£y  -  3Vly/3x,  V'x  -  SV^/Sx,  -  3+/3x  and 
f"  ■  3*$/3x*.  Substituting  Eqs.  (8)  -  (11)  into  Eq  (1),  one  can  show 
that  for  the  ions 


£~  («*  t+  "ix-T-)*  *^r  ♦"'  *(>*‘f^)  ♦" 


U 


k  vr' 


+  (-ik*  —  +  ~  +  1  TT*  S-)  )  ♦'  +  ("kf  (1  +  — )  ♦  7}L  (12) 

7  u  n  “in  7  at  at 


+  2_!zI*I_-2l  V±)  *  -  o 

n  -  n  -  '  T 

at  at 


and  for  the  electrons 


c  k*  n' 
B  («eHtyVey)  n 


(13) 


where  at  •  at  -  kyViy  +  iV'^  and  the  superacrlpta  ('i"('")  Indicate 
first,  second  and  third  derivatives  with  respect  to  x,  respectively* 

tie  assisM  quasineutrality  and  taka  to#  *  dni*  The  following  equation 
is  then  obtained  frost  Eqs.  (12)  and  (13), 


1  +  (l  +  i  Jjl^)  +  (21+  i  J^E  (-ky  +  $-jL  (21) ‘ 


Q  V'  k  V"  k  V' 

-Y~  +  (“ky  (1  +  -?-)  +  j-^L~  +  jl 

y  ay  7  u  at  at 


k  V  Q  k  V"  v 

♦  1  X%~  -1  si'  - 1  Xa-  sl-ia  •> 

3  3=^ 


v,_  k 

fr^-V'V')'  -  0 

yv«y 


(M) 
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V«  ilaplify  Eq.  (14)  by  as  stating  tha  following  ordering  schaae:  J/bt  <  k  , 

2 


Wln/01  <K  1  *  V'  “  V/L  •  V"  “  V/L  *  kyL  <<  °l/vln  »  kyL  <<  V“  where  L 
la  tha  scale  length  of  the  inhomogeneous  plasaa  boundary  layer.  Equation 

(14)  can  now  be  written  as 


iv.  k  V 

in  y  «y 


(irtdv 


■>)  ♦' 


in 


,  <  ,2  VCVB>  -in  kyn*  ,  n- 

>  V  *w  'w  n  n  ~ 

+  M  i.ii.4  \i 


in 


a) 


urt-iv 


in 


k  V _ iv,  _ 

y  sy  in  m 


«tlvln  “ 


-  0 

ey  u 


(15) 


where  «  -  «  -  kyVey(x)  -  w  +  ky(cEx(x)/B) .  Equation  (15)  describes  the 
two-dlaensional  node  structure  of  4  for  the  E  x  B  instability  in  a 
velocity  aheared  plasaa  for  arbitrary  vln/<*> 


III.  RESULTS 


A.  Analytical  Results 

In  general ,  Eq.  (15)  requires  e  numerical  analysis  for  arbitrary 
density  and  electric  field  profiles.  However,  Insight  into  the  nature  of 
the  J  x  instability  can  be  gained  by  first  considering  several  limiting 


cases. 


1.  Locel  Theory 

He  first  reduce  the  differential  equation,  Eq.  (15),  to  an 

albegraic  equation  by  Baking  use  of  local  theory.  That  is,  we  let 

3/ 3te  ♦  ik  ,  and  assuae  k2L2  »  1  and  k2L2  »  1  where  L  -  (n'/n)-1  is  the 
x  x  n  y  n  n 

scale  length  of  the  density  lnhoaogeneity  evaluated  at  x  *  xQ.  For 

A  A 

slapllcity,  we  also  take  E  -  Be  +  E  e  -  constant.  In  this  limit,  Eq. 

~  XX  7  7 

(15)  becomes 


2  2  vin  &  *  t  (c/®) 

lt*X  ♦  k  L  *  -  —  .  m  0 

n  v  n  ~  '*» 

1  (D  U  +  iV. 


Where  »  ■  »  -  k  V _ and  V  „  »  -  cE/B.  Equation  (16)  has  the  solution 

y  ey  ey  x 


.  ',lll  ..  X _  fe  *  ft  Cc/B)  1/2 

* •  - 1 4  i  +  d  +  4^ 


in  n 


Instability  occurs  when 


K  *  •  l 

It  “Ef —  >  °* 

*  n 


The  growth  rates  of  the  Instability  in  the  strong  and  weak  colllsional 
limits  are,  respectively, 


fe  *  l  (c/B) 


;  vln  »  m 


(19) 


J  vln  «  u  (20) 

with  a  real  frequency  u  -  kyV#  in  each  case.  Note  that  Instability  can 
occur  for  Ey  »  0  aa  long  as  ky  *  0  and  *  0  (Eq.  (18);  see  also  Rtskinen 
and  Oisakov  (1982)*  For  k  »  ky,  one  obtains  the  usual  g  x  £  gradient 
drift  instability  growth  rate  (Linson  and  Workman,  1971)  in  Bq.  (19),  and 
the  so-called  high  altitude  limit  (Oesakow  et  al.,  1978)  of  the  E  x  B 
gredlent  drift  instability  in  Eq.  (20). 


2.  Nonlocal  theory 

In  deriving  Iq.  (16)  the  local  approximation  is  used.  That 
is,  the  dispersion  equation  la  solved  baaed  upon  the  plaana  parameters  at 
a  particular  value  of  x,  say  x  xQ;  usually  where  n'/n  is  a  naxlaun  Which 
len4s  to  aexinun  growth.  If  we  now  assume  \  »  EJt(x)  then  a  sheared  E  x  B 
velocity  flow  arlees  Vy  •  V^x)  «  -  cEx(x)/B.  Applying  local  theory  to 
this  situation,  one  might  expect  that  Eq.  (16)  is  still  valid  with  V#y 
evaluated  at  x0,  l.e.,  Vey  -  Vey(xo).  thus,  Eqs.  (19)  and  (20)  follow 
accordingly,  but  the  real  frequency  is  now  given  by  «r  •  ky  V#y(*fl). 
However,  Parkins  and  Doles  (1975)  have  shown,  both  analytically  and  using 
numerical  simulations,  that  this  is  not  the  case.  Mi  do  not  reproduce 
their  detailed  analysis  here,  but  rather,  point  out  the  important  result 
of  their  work. 


Parkins  and  Doles  (1975)  consider  the  strong  collision 
limit  (vln  »  w)  so  that  ion  inertia  terms  can  be  neglected.  Furthermore, 


they  umh  7  •  (a  J)  ■  0  which  leads  to 


n(x)  Ex(x)  ■  «0  *ox  ■  constant 


where  n  ■  n(x  -  -  •)  end  E _ ■  Efx  •  This  is  evident  fron  Eq.  (7) 

O  PX  X 

by  noting  that  B'(x)  ~  E^/I^  end  vlB  »  «  ~  cEy/BI^.  In  this  Halt,  Eq. 
(15)  redness  to 


♦"  +  U  -  — ?-)]  ♦' 


.  If  .  ’S/.f.  8l , . . 

7  s  «  ey  u 


where  2  -  «  -  kyV#y,  V#y-  -(clox/B)(no/n(x))  and  V£y-  “(cE^/BX^/ntx))'. 
Perkins  and  Doles  (1975)  expend  Eq.  (22)  about  x  ■  xQ  where  xQ  is  the 
position  of  naxlsuss  n'/n  by  taking 


n'/o  -  [l-<Jt-«o>2/D2l/l.ii 


»  I  end  kV  »  1,  end  by  asking  several  variable  changes , 
7  ■  7 

they  solve  Iq.  (22)  analytically.  Ike  Important  conclusion  of  their 
theory  is  that  la  that  the  |  x  }  instability  is  stabilised  when 


V%>  2 


Thus,  the  Influence  of  velocity  shear,  i.e.,  an  lnhonogeneous  la  to 
preferentially  stabilise  the  short  wavelength  nodes,  those  with  kyD  »  1. 


B.  NUaerlcal  Basalts 

In  ordsr  to  sol  vs  Eq.  (15)  nunerically,  and  also  to  gain  Insight 
Into  tha  natura  of  the  solutions,  we  transfora  Eq.  (15).  First,  we  note 
that  Eq.  (15)  la  of  tha  fora 

+  P(x)  ♦'  +  q(x)  ♦  -  0  (25) 

where  p(x)  and  q(x)  are  the  coefficients  of  ♦'  and  4,  respectively.  In  Eq. 
(15).  He  let 

♦  -  ♦  exp(-l/2/*p(s)  ds)  (26) 

Substituting  Eq.  (26)  Into  Eq.  (25),  we  find  that  the  transforaed  equation 

Is 

-  Q(*)  ♦  -  0  (27) 

where 


<X*>  -  -q(x)  +  1/2  p'(x)  +  1/4  (p(x))' 


^in 

P(*)  -5-<l  --—TT2-21) 


a  wflv 


In 


-  (£1  -  <21,*,  (l-^2^2-) 


k_V_ 


k  V*  lv 

n'  y  ey  _in  (  JLfZ.  +  JL2I) 


artvln  w 


aflv. 


In 


V. 

u 


7  artv.  IA  «l  l.kk- I  u 


in 


a  n 


arty 


In 


k  V _ lv 


In 


V'  n'  a 


a  n  V-y  n  a 


artv. 


In 
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(28) 


(29) 


(30) 


(31) 
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Equation  (27)  has  a  ala  pie  fora,  albeit  Q(x)  la  a  complicated  function  of 
x,  which  allows  physical  Insight  into  tbs  nature  of  the  node  structure, 
is  an  example,  if  Q  is  real  and  has  Q  >  0  for  |x|  >  xQ  and  Q  <  0  for  |x|  < 
x0  then  one  would  expect  a  bounded  solution  of  +  in  the  region  |x|  <  xQ 
that  exponentially  decays  for  |x|  >  xQ. 

He  now  solve  Eq.  (27)  numerically  for  a  variety  of  conditions  to 
better  understand  tha  influence  of  an  Inhomogeneous  electric  field  on  the 
|><l  instability.  In  all  of  tha  cases  presented,  the  following  density 
profile  is  assumed 

.(»  -  l  ±  (32) 

where  0  <  *  <  1,  L  characterises  the  width  of  the  boundary  layer,  and 
nc  -  n  (x  •  -  •).  By  varying  c,  the  magnitude  of  the  density  gradient 
scale  length  (t#  »  (n rVn)~l)  can  be  changed.  That  Is,  as  c  0, 

I*n  ♦  •  (a  constant  density  profile) ;  as  c  ♦  1,  the  value  of  Ln  ♦  0  (a 
rapidly  changing  denalty  prof 11a).  He  aasune  t  -  0.95  for  the  results 
presented  so  that 

(5l)  -  1.43  L~l  at  x/1  -  -  0.9  (33) 

max 

The  max  In  cm  growth  rate  of  the  instability  is  expected  to  be 

yu  -  l.«  (V0/L)  (34) 

A 

where  V0«  cIq/B  and  we  have  used  Eq.  (19)  assualng  £  •  Eq  ey. 
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Dm  ambient  electric  field  la  chosen  to  be 

£(*)  -  IX(X)  «X  +  By  *y  05) 

where 

E  (x  ■  -*)  »  E  sin  0  e  +  E  cos  6  e  (36) 

^  o  x  o  y 

so  that  6  ■  ten-1  (E^/Ey)  at  X  -  The  Influence  of  the  x  component  of 
the  electric  field  Is  then  studied  by  varying  6  ,  the  angle  between 

'A 

E  and  ey  at  x  -  -  TWo  fonts  of  E^(x)  ere  considered  In  the  analysis : 

Ex(x)  ■  Eq  sin  6  -  constant  (37) 

end 

Ix(x)  -  Eq  sin  0  (no/n(x))  *  constant  (38) 

These  allow  us  to  contrast  the  effects  of  no  velocity  sheer  and  velocity 
sheer  on  the  Instability*  Ms  cosmsnt  that  Eq.  (38)  Is  an  equilibrium 
solution  which  satisfies  V  •  (n  ^)  »  0  la  the  strong  colllolonel  limit 
vla  »  e  (l.e. ,  Eq.  (7)). 

In  Vlg.  3  we  plot  Y  ■  Y/(V  /l)  vs  k.L  for  0  •  0°  end  90°  end  v  •  v/ 

•  y 

(V0/L)  »  1.0  end  100.0*  idiere  Ek  is  chosen  to  be  constant  (Eq.  (37))  and 

V.  •  cE  /B.  A  general  cownsnt  on  ell  of  the  curves  shown  Is  that  Y  Is  an 
o  o 

Increasing  function  k^L  *  but  T  asymptotes  to  e  constant  value  Independent 
2  2 

of  kyL  for  k‘  L  »  1.  This  Is  consistent  with  the  predictions  of  local 
theory.  The  "standard"  case  Is  0  ■  0°,  that  is,  E  •  Eyey.  and  there  Is  no 
component  of  J,  parallel  to  the  density  gradient.  Por  this  case,  two 
values  of  v  are  chosen:  strong  collisions  (v  ■  100.0)  end  week 


collision*  (v  ■  1.0).  As  Is  expected,  the  growth  rate  is  larger  for  the 
larger  value  of  v  in  the  short  wavelength  regime  (kyL  >  1).  Also,  the 
growth  rate  for  v  -  100.0  at  kyL  -  30  la  Y  -  1.39,  and  ia  still 
increasing,  although  slowly,  as  a  function  of  kyL.  This  value  of  y  agrees 
well  with  the  value  obtained  froa  local  theory  (Y  -  1.43  fron  Eq.  (33)). 
The  growth  rate  for  the  weak  collision  case  v  -  1.0  asymptotes  to  a 
somewhat  smaller  value  of  y  (Y  ■  0.79).  However,  note  that  the  difference 
between  the  growth  rates  for  the  strong  and  weak  collislonal  cases  becomes 
smaller  as  kyL  ♦  0,  and  that  the  growth  rates  are,  in  fact,  comparable  for 

A  A 

k  L  *  0.1.  The  "non-standard”  case  is  6  a  90  ,  or  E  ■  E  e  and  the  only 
y  ~  x  x 

component  of  £  is  along  the  density  gradient.  The  major  result  of  this 

limit  is  simply  that  the  instability  can  still  persist  even  though 

Ey  >  0.  The  overall  Influences  of  v  and  kyL  on  the  instability  are  the 

same  as  in  the  previous  case,  8  -  90°. 

In  Fig.  4  ws  plot  Y  vs.  kyL  for  9  »  0°  and  70°  and  v  ■  1.0  and  100.0, 

but  consider  ^  to  be  a  function  of  x  as  in  Eq.  (38)  so  that  velocity 

sheared  flows  occur  for  8  *  0°.  The  curves  for  v  -  1.0  and  100.0 

and  8  a  o°  are  shown  for  comparative  purposes.  The  Important  results  in 

this  figure  are  as  follows.  First,  the  node  is  stable  for  kyL  >  12  for 

both  the  strong  and  weak  collislonal  cases  when  8  -  70° .  This  is  in 

agreement  with  the  conclusion  of  Firkins  and  Doles  (1975);  velocity  shear 

effects  tend  to  stabilise  the  short  wavelength  nodes,  those  such 
2  2 

that  ky  L  »  1.  The  influence  of  shear  on  the  long  wavelength 
nodes  (kyL  <  1)  is  weakly  stabilising.  Second,  the  difference  in  the 
growth  rate  curves  for  v  -  1. 0  and  100.0  is  much  less  than  that  of  the 
case  of  no  shear  (i.e.,  8  ■  0).  And  finally,  since  velocity  shear  can 
stabilise  short  wavelength  modes  before  long  wavelength  nodes,  velocity 


Pig.  4  Plot  of  Y  ■  ‘T/(Vo/L)  ve.  kyL  for  0  ■  0°  and  70°,  and 

for  v  ■  vln/(VQ/L)  *1.0  And  100.0.  The  electric  field  Eg  la 
aesuaed  to  be  Inhoaogeneous  (6q.  (38)). 
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node?  1b  shed  light  on  this  question,  we  consider  the  following 
simplified  equation 

k  (cE  /B)  „ 

♦"  -  [k  -  ik  - \  J  >  T  —  ]  ♦  -  0  (39) 

That  Is,  we  consider  the  llait  v  »  1  and  retain  the  x-dependent ,  Doppler- 

shifted  frequency  (w  -  ky  V#y(x))  as  the  only  contribution  of  the 

inhomogeneous  electric  field  profile*  We  neglect  terms  proportional  to 

V'  ,  v;;,  and  n".  We  emphasise  that  Eq.  (39)  is  not  the  complete  mode 

structure  equation,  but  la  aolved  and  contrasted  to  the  correct  solution 

In  order  to  isolate  a  aingle  effect  of  the  field  inhomogeneity,  vis*,  the 

x  dependent  reeonance  u>  -  k^  V^x).  In  Fig.  7  we  plot  Y  vs*  k^L  for 

6  ■  70°  and  Eg  la  given  by  Eq  •  (38)  *  The  solid  curve  Is  the  solution  to 

Bq.  (27)  for  v  -  100.0,  while  the  dashed  curve  is  the  solution  to  Eq. 

(39).  Although  there  Is  a  small  difference  between  theae  curves 

for  kyL  <  1,  the  Important  point  Is  that  the  mode  is  stabilised  at 

k  L  *  13  in  both  eases .  Thus ,  the  stabilisation  mechanism  Is  related  to 
7 

the  x  dependent  resonance  w  -  ky  V#y(x),  as  opposed  to  velocity  shear 
effects  associated  with  terms  proportional  to  V'y  and  V''.  This  Is  a  key 
result  of  this  analysis. 

We  now  turn  our  attention  to  the  mode  structure  associated  with  the 
£  x  g,  Instability,  and  the  Influence  of  an  Inhomogeneous  slectrlc  field  on 
its  structure.  Figure  8  Is  a  plot  of  the  density  profile  n(x)/n0  (Eq. 
(32)  with  e  »  .095)  and  the  electric  field  profile  Ex(x)/80x  (Eq.  (38) 
with  EQX  «  Eo  sin  8)  versus  x/L.  For  8-0°,  the  electric  field  profile 
la  almply  Eg(x)/^x  -  0.  We  present  plots  of  Q  and  ♦  vs.  x/L  for  these 
profiles.  In  the  subsequent  plots  of  Q  and  ♦,  the  subscript  r  denotes  the 


88 


3 

si 

ft 

(j.i 


SL 


Fig.  7  Plot  of  Y  -  Y/(V0/L)  VI  kyL  for  v  -  vln/(VQ/L)  «  100.0 

and  0  *  70°  using  Bq.  (27)  (solid  curvs)  and  Bq .  (39)  (dashed 


curve).  The  slsctric  field  Eg  is  assuaed  to  be  lnhoaogeneoua 

(Bq.  08)). 


real  part  of  Q  or  ♦,  and  1  denotes  the  Imaginary  part  of  Q  or  +.  the 
paraaeters  considered  for  the  first  set  of  nodes  are  k^L  » 
10.0,  v  -  100.0,  and  9  -  0°  (Fig.  9)  and  6  -  70°  (Fig.  10).  Theae  nodes 
are  considered  to  be  short  wavelength  nodes  since  kyL  »  1. 

Figure  9  is  a  plot  of  Q  (Tig,  9a)  and  9  (Fig.  9b)  versus  x/L  for  the 

case  of  no  electric  field  lnhonogenelty  (8  -  0°  or  Eg(x)  ■  0).  the 

elgenf requency  for  the  node  is  «r  -  0. 0  and  y  •  1. 329.  the  laportant 
points  to  note  are  the  following.  First,  the  wave  potential  Q  is  real  and 
is  such  that  Q  <  0  for  -1.3  <  x/L  <  -0.3  and  Q  >  0  otherwise.  Second,  the 
wave  potential  Q  achieves  a  alnlaiai  at  x/L  *  -0.9,  the  position  of  aaxlata 
1^  (Eq.  (33))*  Third,  consistent  with  this  fora  of  Q,  the  wave  function 
♦  is  a  bounded  node  centered  about  x/L  “  -0. 9  that  falls  off  exponentially 

for  x/L  >  -0.5  and  x/L  <  -1.3.  And  finally,  the  wave  function  la 

reasonably  broad  in  that  its  half -width  at  half  naxlaus  (An)  la  comparable 
to  the  width  of  the  boundary  layer,  i.e..  Ax  «  L/2. 

Figure  10  is  a  plot  of  Q  (Fig.  10a)  and  ♦  (Fig.  10b)  versus  x/L  for 
the  case  of  an  lnhonogeneous  electric  field  (8  -  70°  and  ^(x)  la  shown  in 
Tig.  8).  The  eigenfrequency  for  this  case  la  «r  -  0.5307  and  Y  ■  0.0716. 
Note  that  the  node  has  a  real  frequency  in  contrast  to  the  previous  case 
and  that  the  growth  rate  is  saaller.  Other  laportant  differences  between 
this  situation  and  the  previous  one  are  as  follows.  First,  the  wave 
potential  Q  is  shifted  to  a  larger  value  of  x/L.  The  position  of  the 
ninlaun  value  of  is  at  x/L  “  -0.12.  Also,  note  that  Q  also  has  an 
laaglnary  component .  Second,  the  wavs  function  ♦  is  the  lowest  order  node 
and  has  considerably  aore  structure  in  x/L  than  the  no  shear  case. 
Finally,  the  spatial  extent  of  $  is  sonswhat  narrower  with  Ax  *  0.1  L. 
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A  longer  wavelength  node  Is  now  considered  •  He  choose  k^L  ■  0.1  so 
that  kyL  «  1,  hut  still  consider  v  •  100*0  es  in  the  short  wavelength 
esse.  Figure  11  is  s  plot  of  the  density  profile  n(x) end  electric 
field  profile  Ez(x)/E^)z  for  the  sswe  parameters  es  in  Fig.  8.  However* 
the  range  of  x/L  is  expended  for  conparison  to  the  broadened  node 
structure.  Figure  12  is  a  plot  of  Q  (Fig.  12a)  and  7  (Fig.  12b)  for  the 
esse  of  no  electric  field  inhowogenelty  (•  ■  (P  or  E^(x)  -  0).  Ihe 
elgenfrequency  is  ^  ■  0.0  and  7  *  0.0930.  Ihe  character  of  Q  is 
considerably  different  fron  the  short  wavelength  case  (Fig.  9a).  Ihe 
position  of  the  alnlaun  of  the  potential  well  is  shifted  to  x/L  *  0.0. 
Moreover*  a  “potential  anti-well"  exists  for  -5.0  <  x/L  <0.0  which  tends 
to  inhibit  node  penetration  in  this  region.  Hie  corresponding 
eigenfunction  7  (Fig.  12b)  is  also  substantially  different  fron  the  short 
wavelength  case  (Fig.  9b).  First*  the  wave  function  has  a  surface  wave 
character  la  that  7  *  70  exp  (-«).  Second,  the  wavefunctlon  is 
asyuaetrlcal  about  the  position  of  ninlawa  Qp,  x/L  *  0.0.  Hie 
wavefunctlon  falls  off  very  rapidly  la  the  region  -4.0  <  x/L  <0.0  which 
is  due  to  the  "potential  anti-well"  of  Q  in  this  region.  For  x/L  <  -4.0* 

♦  falls  off  uore  gradually,  slallar  to  its  behavior  for  x/L  >  10.0.  And 
finally*  the  wave  function  is  very  broad*  extending  out  to  x/L  •  50.0. 

Figure  13  is  a  plot  of  Q  (Fig.  13a)  and  7  (Fig.  13b)  versus  x/L  for 
the  eame  parameters  as  Fig.  10*  but  now  we  take  •  ■  70°  so  that  the 
electric  field  is  inhoaogeneous  (see  Fig.  11).  Ihe  elgenfrequency  is 
sr  ■  0.0057  and  7  ■  0.0314.  Although  both  the  wave  potential  Q  and  the 
wave  eigenfunction  7  now  have  laaglaary  components,  Q  and  7  are  quite 
slallar  to  the  no  shear  case.  Ihe  wave  function  is  centered  about  x/L  - 
0.0*  has  an  asynnetrlcal  nature*  and  extends  up  to  x/L  -  50.0.  Thus,  the 
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Mg.  12  Wave  potential  Q  and  «av«  aiganf unction  ♦«*  *  function  of  x/L. 

Tba  subscript*  r  and  i  denote  real  and  imaginary,  reepectively. 
The  parameters  conaldared  are  kyL  -  0.1,  v  •  Vin/(V0/L)  “  100.0 
and  6  •  0^  (i.a.,  ^  -  0).  Hie  eigenfrequency 
la  a  -  0.0  and  Y  -  0.0930.  (a)  Q  vs.  x/L.  (b)  ♦  va.  x/L. 


Fig.  13  Wave  potential  Q  and  vava  eigenfunction  f  a*  a  function  of  x/L 
The  aubacrlpta  r  and  1  denote  real  and  imaginary,  respectively 
The  peraaetere  conaldered  ara  kyL  ■  0. 1,  v  •  vln/(Vo/L)  -  100* 
and  8  -  70°  where  ^  la  given  by  Eq.  (38)*  The  elgenfrequency 
la  »r  «  0*0057  and  7  ■  0.0314.  (a)  Q  vs.  x/L.  (b)  ♦  vs.  x/L. 


Influence  of  the  electric  field  lnhomogeneity  on  the  neve  structure  in  the 

long  wavelength  reglue  (kyL  «  1)  is  much  less  pronounced  then  that  in  the 

short  wavelength  regime  (kyL  »  1).  However,  the  electric  field 

lnhoaogeneity  does  reduce  the  growth  rate  of  the  node  significantly. 

finally  we  present  Fig.  14  which  is  a  marginal  stability  curve 

(i.e.,  Y  »  0)  of  9  vs.  kyL  where  we  have  taken  v  ■  1.0  (dashed  curve) 

end  and  100.0  (solid  curve).  Modes  are  stable  (y  <  0)  and  unstable 

(Y  >  0)  above  and  below  each  of  the  curves,  rsspectlvely.  The  ratio 

Bx(xo)/Ey  is  for  the  case  v  »  100.0  and  has  the  following  weaning.  It  is 

the  ratio  of  E  to  E_  evaluated  at  s  ■  i  ,  Where  r.  Is  the  position  of 
x  y  o  o  w 

node  localisation.  The  position  of  node  localisation  is  defined  as  the 
portion  of  the  nlnlatei  value  of  which  corresponds  to  the  naxiauw  velue 
of  4  .  For  the  parameters  chosen,  it  is  found  that  xQ  -  0.  The  marginal 
stability  criterion  is  then  given  by 


<*o> 


E 


>  0.05  tan  9 


ms 


(40) 


where  9  is  the  value  of  9  at  marginal  stability  and  me  have  used  Eqs. 
(32)  and  (38)  with  e  -  0.95.  From  Fig.  14  we  find  that 


9m  -  88  -  1.4  kyL  (41) 

where  is  measured  in  degrees.  Substituting  Eq.  (41)  into  Eq.  (40),  we 
obtain 


0.05 

tan"a.4kyL) 


(42) 
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Fig.  14  Marginal  atablllty  curva  of  9  vs.  kyL  for  v  -  1.0  (dashad  curva) 
and  v  -  100.0  (aolid  curva).  Tha  aoda  la  atabla  (y  <  0)  and 
unatabla  <Y  >  0)  abova  and  balow  aach  of  thoaa  curvaa, 
respectively  Tha  ratio  BJt(*0)/Ky  va.  kyL  la  for  tha 
caaa  v  -  100.0  uhere  z0/L  “  0.0  ia  tha  poaltlon  of  aoda 
localisation.  In  all  of  thaaa  curvaa,  Eq.  (38)  has  boon  usad  for 


,‘?S‘:V-VV 


for  kyL  «  36*0.  Note  that  Eq.  (43)  is  qualitatively  consistent  with  the 
result  of  Perkins  and  Doles  (1975)  in  that  there  is  an  Inverse 
relationship  between  E  (x  )/E  and  k  .  Also,  Eq.  (43)  is  also 

m.  j  J 

quantitatively  consistent  (see  Eq.  (24))  since  D  >  L  for  the  profiles 
used.  Finally,  as  6  approaches  90°,  i.e.,  Ey  +  0,  the  wavenumber  of  the 
last  unstable  mode  approaches  0.  There  is  no  instability  at  6  -  90° ;  this 
has  been  demonstrated  analytically  by  Perkins  et  al.  (1973). 


IV.  DISCUSSION 


We  have  presented  a  general  theory  of  the  E  x  b  instability  allowing 

for  an  arbitrary  (1)  density  profile,  (2)  inhomogeneous  electric  field 

parallel  to  the  density  gradient,  and  (3)  ratio  of  the  collision  frequency 

to  the  eigenfrequency  (i.e.,  v^/w).  A  differential  equation  is  derived 

which  describes  the  structure  of  the  mode  in  the  direction  of  the 

Inhomogeneity,  which  we  have  considered  to  be  the  x  direction.  The  theory 

is  restricted  to  wave  numbers  such  that  k  L  «  O./v.  and  k  L  «  fl./w: 

y  i  in  y  I 

since  it  has  also  been  assumed  that  v^q/ «  1  and  u/n^«  1  this 
restriction  is  not  Important .  This  work  is  basically  an  extension  of  the 
analysis  of  Perkins  and  Doles  (1975),  whose  theory  is  restricted  to  the 
regime  v^/w  »  1  and  kyL  »  1,  and  considers  a  specific  density  and 

electric  field  profile.  The  principal  results  of  this  study  are  as 
follows . 

1.  For  a  constant  electric  field  profile,  instability  persists  even 

when  Ey  -  0  (0  -  90°).  In  fact,  instability  also  occurs  for 

By  <  0  (9  >  90°)  when  3n/3x  >  0;  this  is  contrary  to  the  simple  one 

A 

dimensional  result  (i.e.,  k  -  k  e  )  which  requires  E  3n/3x  >  0  for 

~  y  y  y 

A 

instability.  Thus,  two-dimensional  mode  structure  (i.e.,  k  -  k^e  + 

A 

kyey)  is  crucial  to  the  instability  (Eq.  (18))  (Linson  and  Workman,  1971). 

2.  For  an  inhomogeneous  electric  field,  an  inhomogeneous  E  x  B 
velocity  occurs  (V  (x)  -  -  cE  (x)/B)  which  has  a  stabilizing  influence  on 

y  * 

the  mode.  Moreover,  the  short  wavelength  modes  (kyL  »  1)  are 

preferentially  stabilized  over  long  wavelength  modes  (kyL  <  1).  This 
result  is  consistent  with  the  work  of  Perkins  and  Doles  (1975). 


a.  The  functional  fora  of  Eg(x)  Is  not  critical  to  stabilisation 
of  the  node*  In  the  absence  of  any  plasna  sources  or  sinks ,  the  Ion 
continuity  equation  gives  the  equllbrlun  relationship  between  the  density 
(n(x))  and  the  electric  field  Eg(x),  as  given  by  Eq.  (7)  (also,  see  Fig. 
8).  Parkins  and  Doles  (1975)  use  this  relationship  in  their  analysis. 
However,  we  have  considered  other  electric  field  profiles  (e.g.,  E^(x)  « 
n(x)  and  Eg(x)  «  tanh  (x)/n(x)).  We  have  found  that  the  instability  is 
still  stabilised  by  the  velocity  Inhoaogenelty,  again  preferentially 
stabilising  the  shorter  wavelength  nodes,  but  that  the  narglnal  stability 
curves  are  different  fron  Fig.  (13). 

b.  The  node  Is  stabilised  because  of  the  x-dependent  resonance 

u  -  kyV#y(x)  In  Eq.  (15).  Tens  proportional  to  JV#y/ix  and  J^/lx2 
are  not  lnportant  for  stabilisation. 

c.  Perkins  at  al.  (1973)  have  shown  analytically  that  the  node 
Is  stable  for  Ej  -  0  and  Ex  given  by  Eq.  (21).  The  numerical  results 
presented  here  are  consistent  with  this  conclusion.  However,  we  add  that 
as  Ey  ♦  0  then  kyL  ♦  0.  This  Is  clear  fron  Fig.  (14)  by  noting  that 

k  L  ♦  0  as  6  ♦  90°. 

y 

d.  In  general,  it  is  found  that  as  vJn  decreases  the  growth  rate 
of  the  node  decreases;  this  Is  expected  fron  linear  theory  (Eqs.  (19)  and 
(20)).  However,  In  the  case  of  an  inhonogeneous  electric  field,  the 
difference  in  growth  rates  between  the  strong  and  weak  collislonal  Units 
considered  Is  net  significant  (see  Fig.  (4)).  Furthemore,  the 
stabilisation  criterion  Is  not  sensitive  to  v^n  (see  Fig.  (14)). 

Ihese  results  are  applicable  to  the  developnent  of  the  E  x  B 
Instability  in  both  barlun  releases  and  the  high  latlt’  *e  F  region 
ionosphere.  First,  the  lnportant  aspects  of  an  Inhonogeneous  electric 


field  on  barium  cloud  strlatlons  has  been  adequately  addresaed  by  Perkins 
and  Doles  (1975).  In  particular,  they  note  that  (1)  the  beck  aide  of  a 
plasma  cloud  nuat  ateepen  aufflclently  ao  that  It  la  alaoat  one- 

A 

diaenalonal  to  allow  the  node  to  grow  (l.e.,  E  *  Ey  ey)  and  (2)  the 
stabilisation  of  the  node  due  to  ^(s)  nay  explain  why  the  aidea  of  a 
plaaea  cloud  do  not  become  unstable.  Furthermore,  f roa  our  studies,  we 
might  hypothesise  that  the  "free ring"  phenomenon  'n  plasma  cloud 
strlatlons  (see  McDonald  et  al.,  1981)  could  be  due  to  shear  stabilisation 
effects,  since  shear  stabilisation  acts  preferentially  on  short  wavelength 
modes,  l.e.,  kyL  »  1,  but  not  on  long  wavelength  modes,  i.e.,  k^L  <  1. 
Second,  the  role  of  an  Inhomogeneous  electric  field  In 
the  g,  x  £  Instability  can  be  very  Important  in  the  structuring  of  plasma 
"blobs"  observed  In  the  high  latitude  F  region.  Experimental  observations 
(Vickrey  et  al.,  1980;  Tsunoda  and  Vickrey,  1982)  Indicate  structuring  In 
both  the  east-west  end  north-south  directions.  Moreover,  small-scale 
structuring  of  the  walls  of  the  "blobs"  have  also  been  observed  and  is 
attributed  to  the  E  x  J  Instability  and/or  the  current  convective 
instability.  The  plasma  configuration  is  not  well-known  but  the 
morphology  of  the  "blobs"  appears  to  be  very  complex.  Not  only  are  there 
Inhomogeneltles  anticipated  In  the  electric  field,  but  there  are  also 
neutral  wind  effects,  field-aligned  pieman  currents,  and  possible  coupling 
effects  between  the  E  and  F  regions.  A  complete  theoretical  treatment 
Incorporating  these  effects  la  beyond  the  scope  of  this  paper.  However, 
the  results  of  this  analysis  strongly  suggest  that  in  order  for 
the  £  x  g  drift  instability  to  be  a  viable  candidate  for  structuring  In 
the  high  latitude  F  region,  then  the  ambient  electric  field  must  be 
orthogonal  or  nearly  orthogonal  to  the  density  gradient.  A  more  complete 


discussion  of  this  problem  will  be  deferred  to  a  later  report  In  which  the 
Influences  of  a  field-aligned  current  and  neutral  wind  are  incorporated 
Into  this  analysis* 
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P.O.  BOX  5400 
AUOQOERQOE,  NM  87115 
OICY  ATTN  DOC  COR  FOR  D.  SHERWOOD 


EG6G,  INC. 

LOS  ALAMOS  DIVISION 
P.O.  BOX  809 
LOS  ALAMOS,  W  85544 
OICY  ATTN  DOC  COM  FOR  J.  BREEDLOVE 

UNIVERSITY  OP  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 
OICY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 
OICY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
OICY  ATTN  DOC  COM  FOR  Ir-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 
P.O.  BOX  1663 
LOS  ALAMOS,  NM  87545 
OICY  ATTN  DOC  CON  FOR  J.  WOLCOTT 
OICY  ATTN  DOC  CON  FOR  R.F.  TASCHEK 
OICY  ATTN  DOC  CON  FOR  E.  JONES 
OICY  ATTN  DOC  COM  FOR  J.  MALIK 
OICY  ATTN  DOC  CON  FOR  R.  JEFFRIES 
OICY  ATTN  DOC  CON  FOR  J.  ZINN  . 

OICY  ATTN  DOC  COM  FOR  P.  KEATON 
OICY  ATTN  DOC  CON  FOR  D.  WESTERVELT 

SAMDIA  LABORATORIES 
P.O.  BOK  5800 
AIBUQUERQUE,  1M  87115 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DQC  CON  FOR  A.  THORNS ROUGH 

OICY  ATTN  DOC  CON  FOE  T.  WRIGHT 

OICY  ATTN  DOC  COM  FOR  D.  DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 

SAMDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOK  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  MtytFHKY 
OICY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OP  MILITARY  APPLICATION 
DEPARTMENT  OP  ENERGY 
WASHINGTON,  D.C.  20545 
OICY  ATTN  DOC  CON  DR.  TO  SONG 


DEPARTMENT  OP  CONEUCE 
NATIONAL  BUREAU  OP  STANDARDS 
WASHINGTON,  R.C.  20234 

(ALL  CORRESl  ATTN  SEC  OFFICER  FOR) 
OICY  ATTN  R.  MOORE 


INSTITUTE  FOE  TELECOM  SCIENCES 
RATIONAL  TELECOMMUNICATIONS  A  INFO  AOMIN 
BOULDER,  CO  80303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

01CY  ATTN  «.  UTLAUT 
01CT  ATTN  0.  CSOMBIE 
01CY  ATTN  L.  BESET 

NATIONAL  OCEANIC  A  ATMOSPHEEIC  AOMIN 
ENVIEOIMENTAL  RESEAECH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
01CY  ATTN  R.  GRUBB 
01CT  ATTN  AERONOMT  LAB  G.  REID 

DEPARTMENT  OF  DHFEM8E  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
01CY  ATTN  I.  OARFUNEEL 
01CT  ATTN  T.  8AUH 
01CT  ATTN  9.  J08EPHS0N 
OICT  ATTN  S.  BONER 
01CY  ATTN  N.  STOCKHELL 
Old  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  COW 
5  OLD  OONOOID  MAD 
BURLINGTON,  HA  01803 
Old  ATTN  RADIO  SCIENCES 

SEREELET  RESEARCH  ASSOCIATES,  INC. 

P.O.  BOX  983 
BERKELEY,  CA  94701 
Old  ATTN  J.  NOREMAN 
Old  ATTN  C.  PRETTIE 

BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  HA  98124 
Old  ATTN  6.  KEISTER 
Old  ATTN  D.  MURRAY 
Old  ATTN  0.  HALL 
Old  ATTN  J.  KENNEY 

BROWN  ENGINEERING  COMPANY,  INC. 

CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  35807 
Old  ATTN  ROMM  A.  DELIS  ER IS 

CALIFORNIA  AT  SAM  DIEGO,  UNI V  OF 

P.O.  BOX  4049 

SAN  DIEGO,  a  92106 


CHARLES  STARK  DRAPER  LABORATORY,  INC. 

555  TECHNOLOGY  SQUARE 
CAM RIDGE,  MA  02139 
Old  ATTN  D.B.  COX 
Old  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTKICUK  ROAD 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 
Old  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
Old  ATTN  H.  LOGSTOM 
Old  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL,  IRC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
Old  ATTN  J.  ROBERTS 
Old  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
QODGASD  BLVD  KIM  OP  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
Old  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
Old  ATTN  F.  REIBERT 

GENERAL  ELECTRIC  TECH  SERVICES  CO. ,  INC. 
HME8 

COURT  STREET 
SYRACUSE,  NY  13201 
Old  ATTN  G.  HILLMAN 

GENERAL  research  corporation 
SANTA  BARBARA  DIVISION 
P.O.  BOX  6770 
SANTA  BARBARA,  CA  93111 
Old  ATTN  JOHN  ISE,  JR. 

Old  ATTN  JOEL  CARBARINO 


GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01 CY  ATTN  T.N.  DAVIS  (UN CLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
01  CY  ATTN  NEAL  BROW  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01 CY  ATTN  MARSHALL  CROSS 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 
01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 
01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 
01CY  ATTN  J.M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01 CY  ATTN  HANS  NOLFARD 
01CY  ATTN  JOEL  BENGSTQN 

INTL  TEL  A  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 


KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
01CY  ATTN  DASIAC 
01CY  ATTN  TIM  STEPHANS 
01CT  ATTN  WARREN  S.  KNAPP 
01CY  ATTN  WILLIAM  MCNAMARA 
01CY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
01  CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  6  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
01CY  ATTN  DEPT  60-12 
01CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  A  SPACE  CO. ,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 
01CY  ATTN  MARTIN  HALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
01CY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

LOCKHEED  MISSILE  &  SPACE  CO. ,  INC. 
HUNTSVILLE  RESEARCH  &  ENGS..  CTR. 

4800  BRADFORD  DRIVE 
HUNTSVILLE,  AL  35807 
ATTN  DALE  H.  DIVIS 


JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO.  CA  92138 
01  CY  ATTN  J.L.  SPERLING 

JOWS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOWS  HOPKINS  ROAD 
i  ahmi  ,  ib  »707 
OICY  ATM  DOCUMENT  LIBRARIAN 
01  CY  ATTN  THOMAS  POTEKRA 
OICY  ATTN  JOW  OASSOULAS 
OICY  ATTN  DR.  DONALD  J.  WILLIAMS 


MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
OBLANDO,  PL  32805 
OICY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 
OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  P.  WALDRON 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 


MCDONHEL  DOUGLAS  CORPORATION 
$301  BCLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CY  ATTN  N.  HARRIS 
01CT  ATTN  J.  HOULE 
01CT  ATTN  GEORGE  MROZ 
01CT  ATTN  W.  OLSON 
Old  ATTN  R.W.  RALPRIN 
Old  ATTN  TEOMICAL  LIBRARY  SERVICES 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01CY  ATTN  P.  FISCHER 
01CT  ATTN  W.F.  CREVIER 
01CT  ATTN  STEVEN  L.  GUTSCHE 
01CT  ATTN  D.  SAPPENFIELD 
01CY  ATTN  R.  BOGUSCH 
01CT  ATTN  R.  HENDRICK 
Old  ATTN  RALPH  KILB 
Old  ATTN  DAVE  SOULE 
Old  ATTN  P.  FAJHI 
01CT  ATTN  H.  SCHEIBE 
Old  ATTN  CONRAD  L.  LCNGMIRE 
Old  ATTN  WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 
Old  ATTN  JOOt  MORGAN  STERN 
Old  ATTN  G.  HARDING 
Old  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
Old  ATTN  W.  HALL 
Old  ATTN  W.  POSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CL OVERFIELD  BLVD. 

SANTA  MONICA,  CA  90404 
Old  ATTN  B.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  a ASS  TO  THIS  ADDRESS) 

Old  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 
Old  ATTN  IRVING  L.  KOFSKY 


PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 
Old  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
Old  ATTN  FORREST  GILMORE 
Old  ATTN  BRYAN  GABBARD 
Old  ATTN  WILLIAM  B.  WRIGHT,  JR. 
Old  ATTN  ROBERT  P.  LELSVIBR 
Old  ATTN  WILLIAM  J.  KARZAS 
Old  ATTN  H.  ORY 
Old  ATTN  C.  MACDONALD 
Old  ATTN  R.  TURCO 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
Old  ATTN  CULLEN  CRAIN 
Old  ATTN  ED  BEDIOZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 
Old  AXIS  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  BID  AVENUE 
NEW  YORK,  NY  10023 

Old  ATTN  VOICE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.O.  BOX  2351 
LA  JOLLA,  CA  92038 
Old  ATTN  LEWIS  M.  LIN  SON 
Old  ATTN  DANIEL  A.  HAMLIN 
Old  ATTN  E.  FRIEMAN 
Old  ATM  E.A.  STRAKER 
Old  ATTN  CURTIS  A.  SMITH 
Old  ATTN  JACK  HCPOMCALI. 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDCE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 
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SRI  INTERNATIONAL 
333  RAVENS WOOD  AVENUE 
MENLO  PARK,  CA  94025 
01  CY  ATTN  DONALD  NEILSON 
OICY  ATTN  ALAN  BURNS 
01CT  ATTN  G.  SMITH 
OICY  ATTN  L.L.  COBB 
OICY  ATTN  DAVID  A.  JOOiSON 
OICY  ATTN  WALTER  G.  CHE SHUT 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  M.  BARON 
OICY  ATTN  RAY  L.  LEADABRAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  J.  PETERSON 
OICY  ATTN  R.  HAKE,  JR. 

OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 

STEWART  RADIANCE  LABORATORY 
UTAH  STATE  UNIVERSITY 
1  DE  ANGELO  DRIVE 
BEDFORD,  MA  01730 
OICY  ATTN  J.  ULWICK 

TECHNOLOGY  INTERNATIONAL  CORF 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.P.  BOQUIST 

TRW  DEFENSE  A  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 

VISIDYNB 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPBfTER 
OICY  ATTN  C.  HUMPHREY 


IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 


OF  THE  FOLLOWING  PEOPLE 


PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 
DR.  P.  MANGE  -  CODE  4101 
DR.  R.  MEIER  -  CODE  4141 
DR.  E.  SZUSZCZEWICZ  -  CODE  4187 
DR.  J.  GOOEMAN  -  CODE  4180 
DR.  R.  RODRIGUEZ  -  CODE  4187 
CODE  2628  -  20CY 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANS  COM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  J.  AARONS 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW1 6-262 
CAMBRIDGE,  MA  02139 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
DR.  S.  CHANDRA 
DR.  K.  MAEDA 
DR.  R.F.  BENSON 

NATIONAL  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 
12CY  ATTN  TC 

COMMANDER 

NAVAL  AIR  SYSTEMS  OOtttAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO.  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL 
LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 
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PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08540 
DR.  F.  PERKINS 

SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  L.  LINSOM 
OR.  E.  FRIEMAM 

STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  RANKS 


UNIVERSITY  OP  PITTSBURGH 
PITTSBURGH,  PA  15213 
OR.  N.  ZABUSKY 
OR.  M.  BXONDI 
OR.  E.  OVERMAN 

VTAN  STATE  UNIVERSITY 
4TR  AND  OTH  STREETS 
1RRAR,  STAN  S4322 
OR.  R.  HARRIS 
OR.  R.  BAKER 
NL  R.  SCMUMK 


y 


UNIVERSITY  OF  CALIFORNIA, 
LOS  ANGELES 
405  HILLGARD  AVENUE 
LOS  ANGELES,  CA  90024 
DR.  F.V.  GOR0NXTX 
DR.  C.  KENNEL 
DR.  A.Y.  WONG 

UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 
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